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Abstract 
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We use both Abel's lemma on summation by parts and Zeilberger's algorithm 
to find recurrence relations for definite summations. The role of Abel's lemma can 
be extended to the case of linear difference operators with polynomial coefficients. 
\ This approach can be used to verify and discover identities involving harmonic 

I numbers and derangement numbers. As examples, we use the Abel-Zeilberger al- 

gorithm to prove the Paule-Schneider identities, the Apery-Schmidt-Strehl identity, 
Calkin's identity and some identities involving Fibonacci numbers. 
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1. Introduction 



^ ■ The main idea of this paper is to employ the classical lemma of Abel and Zeilberger's 

. algorithm for hypergeometric sums to verify and to discover various identities on sum- 

mations that are not hypergeometric. Abel's lemma [T] on summation by parts is stated 
as follows. 

Lemma 1.1 For two arbitrary sequences {ak} and {hk}, we have 

n— 1 n— 1 

^(ofc+i - ak)hk = ^ ak+i{bk - bk+i) + a„6n - am&m- 

k=m k=in 

For a sequence {rk}, define the forward difference operator A by 
Then Abel's lemma may be written as 

n— 1 n—l 

fefcAofc = - ^ Uk+iAbk + anbn - ambm- (1-1) 

k=m k=m 
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Graham, Knuth and Patashnik [TT] reformulated Abel's lemma in terms of finite 
calculus to evaluate indefinite sums. Recently, Chu [8] utilized Abie's lemma to prove 
basic hypergeometric identities including Bailey's very well poised g^e-series identity by 
finding pairs (0^,6^). Applying Abel's lemma to the pairs (0^,6^), one obtains contigu- 
ous relations for the basic hypergeometric sums. Chen, Chen and Gu [B] presented a 
systematic approach to finding pairs (a^, b^) by using the g-Gosper algorithm. 

In this paper, we combine Abel's lemma and Zeilberger's algorithm to verify and dis- 
cover identities. Moreover, we use an extended version of Abel's lemma to deal with sums 
involving holonomic sequences. Let us illustrate our approach by considering identities 
involving harmonic numbers. The k-th harmonic number is given by 



k 



Note that by definition, Hj. = whenever k < 0. Let fk be a hypergeometric term, i.e., 
fk+i/fk is a rational function of k. We focus on the summation 



(1.3) 



k=m 

We can use Gosper's algorithm [TU] to determine whether there exists a hypergeometric 
term such that Aa^ = fk- If such ak exists, by Abel's lemma we get 

n—1 n—1 n—1 

fkHk = ^ HkAak = - ^ + ttnHn - amHm- (1.4) 

k=m k=m k=m 

Hence we can transform a summation involving harmonic numbers into a hypergeometric 
summation. We call such an approach the Abel-Gosper method. 



The same idea applies to the definite summation 

n 

5(n) = ^F(n, k)Hk, 

k=0 

where F{n, k) is a proper hypergeometric term such that -F(?t,, /c) = for k > n. In 
this case, we can apply Zeilberger's algorithm to find a hypergeometric term G(n, k) and 
polynomials ^^^{n), . . . ,Pd{n) such that 

d 

'^Pj{n)F{n + j, k) = G{n, k + 1) - G{n, k). 

j=0 

Taking 

d 

fk = ^Pj{n)F{n + j, k) and = G(n, k) 

j=0 

in (ll.4p and summing over k from to n + rf, we deduce that 

d n+d L I 1 N 

Y,PM)S{n + j) = - E k+l + + + ^)Hn+d+i. (1.5) 

j=0 k=0 
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Moreover, we see that the summation on the right hand side of (11.51) is again a hyper- 
geometric sum. This approach is called the Abel-Zeilberger method. 

In order to apply the above approach to general holonomic sequences, we extend 
Abel's lemma by replacing the operator A with a linear difference operator L of the 
form 

Lak = ro{k)ak + ri{k)ak+i H h rd{k)ak+d, 

where each ri{k) is a rational function of k. By combining a variation of Zeilberger's 
algorithm and the extended Abel's lemma, we can find a recurrence relation for a sum 
of the form 

Sin) = J2fi^^k)gin,k), (1.6) 

k 

where f{n, k) is hypergeometric and g{n, k) satisfies two recurrence relations 
g{n, k + d) = ri(n, k)g{n, k) + ■ ■ ■ + rd{n, k)g{n, k + d — 1) + u{n, k), 

and 

g{n + l,k) = s{n, k)g{n, k) + v{n, k), 

where all the coefficients ri{n,k) and s{n,k) are rational functions and u{n, k),v{n, k) 
are hypergeometric terms. The algorithm for finding recurrence relations for the sum 
(II. 6p is called the Abel-Zeilberger algorithm. 

The paper is organized as follows. In Section 2, we give examples to demonstrate 
that many indefinite sums involving harmonic numbers can be reduced to hypergeomet- 
ric sums by the Abel-Gosper method. Section 3 shows how to apply the Abel-Zeilberger 
method to find recurrence relations of definite sums involving harmonic numbers. For 
example, the Paule-Schneider identities fall into this framework. In Section 4, by extend- 
ing Abel's lemma and using a variation of Zeilberger's algorithm, we present the Abel- 
Zeilberger algorithm. The last section provides several examples of the Abel-Zeilberger 
algorithm including identities involving Fibonacci numbers and derangement numbers, 
as well some identities on multiple sums. 

2. The Abel-Gosper method 

In this section, we give several examples to illustrate how to combine Abel's lemma with 
Gosper's algorithm to evaluate indefinite summations. We shall focus on summations 
involving harmonic numbers. 

We begin with a simple example. Consider the sum 

n 

S{n) = J2Hk. 

k=l 

Expressing 1 as Ak, we obtain 

n n 

S{n) = J2 HkAk = - + l)AHk + {n + l)Hn+i -Hi = {n + l)Hn - n. (2.1) 

k=l k=l 
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The idea of the above example can be generahzed to indefinite sums of products of 
polynomials and harmonic numbers. The following result is due to SpieB [15]. Here we 
give a derivation based on Abel's lemma. 



Theorem 2.1 Let Hk be the k-th harmonic number and u{n) be a polynomial of degree 
m in n. Then 

n 

J2u{k)Hk=p{n)Rn-q{n), n = 0,1,2,..., (2.2) 

fc=0 

where p{n) andq{n) are both polynomials of degree m+1 inn. Moreover, p{n) is divisible 
by n + 1. 

Proof. It is well-known that there exists a polynomial f{k) of degree m + 1 in k such 
that Af{k) = u{k) and the constant term of f{k) is zero. Therefore, we may write 
f{k) = kg{k), where g{k) is a polynomial of degree m. From (11.41) it follows that 

n n n 

<k)Hk = J2 ^{kg{k))Hk = - ^ (^(A; + 1) + (n + l)g{n + l)i7„+i. 

fc=0 fc=0 fc=0 

Since the sum X]fc=ofi'(^ + 1) is a polynomial of degree m + 1 in n and 
(n + l)g{n + l)i/„+i = {n + l)g{n + l)Hr, + g{n + 1), 
we arrive at (12. 2p . I 
Setting u{n) = 1 in (12.21) . we obtain (12. ip . When u{n) = for m = 1, 2, 3, we have 

^,ft = !!(!!±i)ff„-(!i_l)!!, (2.3) 

k=l 



n 



y _ _ ("-lM-"' + l). (2.4) 

^-^ 6 36 

fc=i 

fc=l 

The same idea also applies to the bonus problem 69 proposed by Graham, Knuth 
and Patashnik [HI Chapter 6]. 

Example 2.2 Find a closed form for 



n 

' 1.2 TJ 

n+k- 



k=l 

The above sum can be rewritten as 

2n 2n 



J] (A; - nfHk = ^(A; - nfHk - Y,{k - ufH^. (2.6) 



k=n+l k=l k=l 
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Expanding the summands and applying formulas (12. ip . fl2.3p . (12 ■4p . we obtain that 

± Pff„,, . "("+i)p"+i) pg,„ _ _ (2.7) 

fe=l 

We may also apply Gosper's algorithm directly to {k — n)^ and get its indefinite sum 

1 , 3 272 + 1,2 6n2 + 6n + l, ^ 

-F k^ H A; + C. 

3 2 6 

Setting C = and using fll.4p . we also obtain (12. 7p . 

We remark that Chyzak [9] and Schneider [T3l[Tl] have proved (12. 7p by an extension 
of Zeilberger's algorithm and Karr's algorithm, respectively. 



Example 2.3 Evaluate the sum 

k=0 



1 f2k 



^2k f2k^ _ 1 (2k 



By Gosper's algorithm, we find that 

2k /2 

'^\k J A^\kJ' 

Therefore, 

n—l ^ /„ , \ n— 1 



k — k — 



^ n + 1 /2n + 2\ 2n/2n\^ 

Example 2.4 We have 

n 

Hi = {n + l)Hl - (2n + l)Hn + 2n. {21 

k=0 

Proof. Setting = k and bk = in (11. ip . we deduce that 

n n 

5^ = - Y,{k + 1) Ai/^ + + 



k=0 k=0 



n n ^ 

k=0 k=0 

-2{n + 1)H^ + 2n- + (n + l)i7^^i (By (O) 

(n + l)iJ,^ _ (2n + l)if„ + 2n. I 
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Similarly, by setting bk = H^, we deduce the following identities, see [9|[T5]. 

n 

J2 Hi = {n + l)Hl - 2 (2^ + + 3(2n + + -H^^^^ - 6n, (2.9) 

k=0 

Y,{2k + l)Hl = {n + If Hi - ln{n + l)Hl + g„ - + 1). (2.10) 

fc=0 

3. The Abel-Zeilberger method 

In this section, we use Abel's lemma and Zeilberger's algorithm to find recurrence rela- 
tions of the definite summation of the form 

n 

S{n) = Y,F{n, k)Hk, 

k=0 

where F{n, k) is a hypergeometric term in two variables n and k. We shall give two 
examples to explain the method. First, we consider an identity due to Chu and De 
Donno [7j 

Example 3.1 For n > 0, we have 

x:QV.=(2i/„ -//,.)(';)• (3.1) 

Proof. By applying Zeilberger's algorithm to (^)^, we obtain a skew recurrence relation 
+ "^^^'' + ^^(1) =G{n,k + l)-G{n,k), 

where 

G{n, k) = {-3-3n + 2k) 

Let 

Substituting F{n,k) = (^Y and G{n,k) = (— 3 — 3n + 2/;;) into (11. 5p . we find that 

(„ + l)S(„ + l)-2(2n+l)S(n)=g^^(J =— (J. (3^2) 

where the second equality can be justified by applying Zeilberger's algorithm. 
It is easy to verify that 

R{n) = {2Hn - H2n) 

6 



G-0 




satisfies the same recurrence relation ( 13.21) . Since S{0) = -R(O) = 0, we get (13.11) . This 
completes the proof. I 



In [12], Paule and Schneider considered the following summations: 

a = 1,2,.... (3.3) 



fc=0 

They found closed forms of T'^f^ for 1 < a < 4 and derived recurrence relations of T^"^ 
for 5 < a < 9. As will be seen, we can combine Abel's lemma and Zeilberger's algorithm 
to deal with the summations . As an example, let us consider the case a = 3. 

Example 3.2 For n > 0, we have 

) = (-l)^ 

Proof. Let 

F{n, k) = {n- 2k) I 
By Zeilberger's algorithm, we find that 

F(n, k) + F{n + l,k) = G{n, k + 1) - G{n, k), 

where 

G{n, k) = {2n-k + 2) 

Let 



, 3 

n 



, 3 

n 



k-l 



S{n) = Tf = + 3(n - 2k)H,) Q . 



fe=0 

By (II. Sp . we deduce that 

.(„).5(..i)^-3E^(:)%E(:)%ECtT- 

k=0 ^ ^ k=0 ^ ^ k=0 ^ ^ 

By Zeilberger's algorithm, we find that the right hand side, denoted by R{n), satisfies 

(n + l)R{n) + (n + 2)R{n + 1) = 0. 

Since -R(O) = 0, we have R{n) = for n = 0, 1, . . .. It is clear that S{0) = 1. So we get 
S{n) = (—1)". This completes the proof. I 

Moreover, as a direct consequence of (II. Sp . we have the following property. 
Theorem 3.3 Let 

f^'=E("-2*)(; 
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Assume that the minimal recurrence relation for ulf^ computed by Zeilberger's algorithm 
is 

1=0 

Then the summation 

d 

T.P^in)Tit (3.4) 

is a hypergeometric summation. 

We find that the sum (13 ■4p equals zero for a = 1, 2, . . . , 9 and conjecture that it holds 
for any nonnegative integer a. We note that this conjecture implies the conjecture of 
Schneider and Paule [12], which says that T^'^^ satisfies the minimal recurrence relation 
for Un"'^ computed by Zeilberger's algorithm. 

4. The Abel-Zeilberger algorithm 

In this section, we give a description of the Abel-Zeilberger algorithm. Notice that in 
the applications of Abel's lemma given in previous sections, the main idea lies in the fact 
that AHk is a hypergeometric term. In fact, there are other sequences satisfying similar 
properties that lead us to consider an extension of Abel's lemma. 

Let {flfc} be an arbitrary sequence. We consider a linear operator L of the form 
Lofc = ro{k)ak + ri{k)ak+i H h rd{k)ak+d, 

where each rj{k) is a rational function of k. We associate the operator L with a dual 
operator L* defined by 

L*ak = ro{k)ak + ri{k - l)afc_i H h r^i^k - d)ak-d- 

In the above notation, Abel's lemma can be extended to the following form. 

Lemma 4.1 For two arbitrary sequences {ak} and {bk}, we have 

n—l n— 1 

L*ak ■bk = Y,^k- Lbk - T{n) + T(m), (4.1) 

k=m k=m 

where 

d i 

= ^^r,(/c-j>fc_,6fc+,_,-. (4.2) 
i=i j=i 
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Proof. It is easy to verify that 



n-l 



n-1 



k=m 



^ L*ak ■ hk = "^"^rii^k - i)ak-ihk 

k=m i=0 
d n—l—i 

= ^ ^ ri{k)akhk+i 



j=0 k=m- 
n-1 



^ro{k)akbk + ^ 



k=m 
n-l 



1=1 



n— 1 m — 1 n— 1 



k=m k=m—i k=n—i. 



ri{k)akbk^ 



d m—1 



d n—1 



^ afc ■ + ^ ^ ri{k)akbk+i - ^ ^ ri{k)akbk+i 



k=m 



i=l k=m—i 



1=1 k=n- 



Let f{n,k) be a bivariate hypergeometric term and g{n,k) be a bivariate function. 
We aim to find a linear recurrence relation for the definite sum 



k=m 

Suppose that there exist rational functions rj{n,k) such that 

d 

Lg{n, k) = ^rj(n, k)g{n,k + j) 

j=0 

is a bivariate hypergeometric term. We shall try to find polynomials Pi{n), which are 
independent of k and not all zero, together with hypergeometric terms a{n, k) and w{n., k) 
such that 



^^Pj(?T.)/(n + i, k)g{n + i,k) = g{n, k)L*a{n, k) + w{n, k). (4.3) 
Summing f l4.3p over k and applying the extended Abel's lemma, we deduce that 

^Pi{n)S{n + = X] X] Pi('^)f(^ + ^' + i, k) 

i=0 k=m 

= ^^((^(n, k)L*a{n, k) + /c)) 

k=m 

= a{n, k)Lg{n, k) + Y^ w{n, k) - T{i + 1) + T(m), (4.4) 



i=0 



k=m 



k=m 



where T{k) is given by (14. 2p . Notice that in the last expression, the two summands 
are hypergeometric. By Zeilberger's algorithm, the two sums satisfy linear recurrence 
relations, which lead to a non-homogenous linear recurrence relation for S{n). 
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To solve equation (14. 3p . we impose the condition that g{n, k) satisfies the relation 

g{n + l,k) = s{n, k)g{n, k) + v{n, k), 

where s{n, k) is a rational function and v{n, k) is a hypergeometric term. By induction, 
it is easy to show that there exist rational functions Si{n,k) and hypergeometric terms 
Vi{n, k) such that 

g(n + i,k) = Si{n, k)g{n, k) + Vi{n, k). (4.5) 

Now we can solve the following equation for Pi{n) and a{n,k) by a variation of 
Zeilberger's algorithm 

I 

^^Pi(n)/(?7, + i, k)si{n, k) = L*a{n, k). (4.6) 

It can be seen that a(n, k) is similar to /(n, /c), that is, 

a(n,k) 
f{n,k) 

is a rational function of n and k. Hence (14.61) is equivalent to 

Since f{n,k) is hypergeometric, both f{n + i,k)/f{n,k) and f{n,k — j)/f{n,k) are 
rational functions. Therefore, (14. 7p is a non-homogenous linear recurrence equation on 
R{n, k) with parameters Pi{n), which can be solved by Abramov's algorithm [2]. 

Once we find a solution {po{n), . . . ,pj{n), a{n, k)) to (14. 6p . it is easy to check that 

{Poiji), . . .,pi{n), a(n, k),w{n, k)) 
is a solution to (14.30 . where 



i=0 



k) = J2Pi(^)fi^ + ^^^)^^i^^^)- (4-8) 



In summary, the Abel-Zeilberger algorithm can be described as follows. 

Input: a hypergeometric term f{n,k) and a term g{n,k) satisfying two recurrence 
relations 

g{n, k + d) = ri(n, k)g{n, k) + ■ ■ ■ + rd{n, k)g{n, k + d — 1) + u{n, k), (4.9) 

and 

g{n + l,k) = s(n, k)g{n, k) + v{n, k), (4-10) 

where ri{n, k) and s{n, k) are rational functions, and u{n, k) and v{n, k) are hyper- 
geometric terms. 
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Output: polynomials Pi{n) that are independent of k, two hypergeometric terms ti{n, k), 
t2{n, k) and a term T{k) satisfying 

^Pi{n)S{n + i) = ^ti(n,A;) + ^t2(ri,A;)-T(£+l) + r(m), 

i=0 k=m k=m 

where 

^^'^) = X] f('^^^)9in,k). 

k=m 

The algorithm consists of the following steps. 
Initially, we set / = 0. 

Step 1. For < i < I, compute the rational functions Si{n,k) and the hypergeometric 
terms Vi{n,k) defined by fl4.5p by using the recurrence relations 

Si+i{n, k) = s{n + i, k)si{n, k), 

Vi+i{n, k) = s{n + i, k)vi{n, k) + v{n + i,k), 

with the initial values so{n, k) = 1 and vo{n, k) = 0. 
Step 2. Let 

Lg{n, k) = — ri(n, k)g{n, k) — ■ ■ ■ — r^in, k)g{n, k + d — 1) + g{n, k + d). (4.11) 

According to f l4.7p . construct an equation on Piin) and R{n, k) = a{n, k)/ f{n, k). That 
is, compute polynomials 

Pj{n,k),0<j <d and Qiin,k), < i < I 

such that 

d I 

Pjin, k)R{n, k-j) = 5^p.(n)g,(n, k). (4.12) 

Step 3. Solve equation (14.12p for R{n,k) and pi{n),0 < i < I hj using Abramov's 
algorithm. If all the polynomials Pi{n) are zeros, then we increase / by one and repeat 
steps 1-3. 

Step 4. Compute w{n, k) based on 04.81) . For L given by f l4.1ip . compute T{k) according 
to do}. Finally, set 

ti{n, k) = a{n, k)u{n, k) and t2{n,k) = w{n,k). 

Then {pi{n),ti{n, k),t2{n, k),T{k)) is the desired output. 
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5. Examples 



In this section, we provide several examples to compute summations by using the Abel- 
Zeilberger algorithm. 

We first consider the case when g{n, k) is independent of n in the Abel-Zeilberger 
algorithm as described in the previous section, i.e., s(n, k) = 1 and v{n, A;) = in (14.101) . 
In this case, we have Si{n, k) = 1 and Vi{n, k) = so that t2{n, k) = 0. 

Let Fn be the n-th Fibonacci number which is defined by the recurrence relation 

Fn+2 = Fn+l + Fn, 

for 77, > 0, with initial values Fq = 0,Fi = 1. Employing the Abel-Zeilberger algorithm, 
we can prove the following identities on the product of binomial coefficients and Fibonacci 
numbers, see [^[T7]. 



Example 5.1 We have 



k=0 



(5.1) 



E(-i)'(:) 



Fk = -Fn, (5.2) 



fc=0 

n 



k=0 

n 



k=0 



^(^)F4fc = 3"F2„. (5.4) 



r F3, = 2"F2„, (5.3) 



Proof. For equation fl5.1l) . taking f{n,k) = (^) and g{n,k) = F^ as the input of the 
Abel-Zeilberger algorithm, we obtain 

Po{n) = 1, pi{n) = -3, p2{n) = 1, ti{n, k) = t2{n, k) = 0, 

and 

Tik) = {i2k - n - m + (n + 2 - k)F,,,)^—^^^^. 

We see that T{n -|- 3) = T(0) = 0. Therefore, the summation 



fc=0 

satisfies the recurrence relation 



S{n) - 3S{n + 1) + S{n + 2) = 0. 

Since F2n satisfies the same recurrence relation with initial values Fq = and F2 = 1, 
we deduce that S{n) = F2„. 
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The other three identities can be proved in the same fashion. The detailed arguments 
are omitted. This completes the proof. I 

It is not difficult to see that the Abel-Zeilberger algorithm could be used to ver- 
ify identities involving general C-finite sequences. For instance, suppose that {Gk}k>o 
satisfies a recurrence relation 

Gk+2 = bGk+i + cGfc, A; > 0, 

where b and c are constants. Then the Abel-Zeilberger algorithm generates the recurrence 
relation 

(6 + 1 - c)S{n) - (6 + 2)S{n + 1) + S(n + 2) = 

for the summation 



i — n V / 



fc=0 

The next example involves the n-th derangement number Dn as given by 

-1^ 



k=0 



k\ 



Using the method of MacMahon's partition analysis, Andrews and Paule [3] the following 
identity on D„. We shall give a derivation by applying the Abel-Zeilberger algorithm. 

Example 5.2 

for n > N > n — k. 

Proof. Substituting k + N — j for j, the left hand side of (15.51) can be rewritten as 

k+N 



j=N 

Because of the recurrence relation 

Dn = uDn^i + (-1)" 

for Dn, we may put 

as the input of the Abel-Zeilberger algorithm. Then we obtain 

13 



Denote the right hand side by G{k). By Zeilberger's algorithm, we find that for A; > 0, 

{l + N -n + k)G{k) -{2 + k + N)G{k + 1) = 0, 

which imphes that G{k) = for k > n — N > 0. Thus we get ^(A^, n) = S{N + 1, n). In 
particular, 



SiN,n) = Sin,n) = J2(^^^^'jD,. 



3>n 

Applying the Abel-Zeilberger algorithm to 



/(ra, j) = ( . ^ ) and g{n,i) = Dj, 
J — n 



we find 



(n + l)S(n, n) + (n + k + l)S(n + l,n + 1) - S(n + 2,n + 2) = Vf-lV 1^1= 0. 

j>n ' 

By Zeilberger's algorithm, we see that the right hand side of (15. 5 p satisfies the same 
recurrence relation. Finally, from the identity 

i=o ^•'^ 

we deduce that S'(0, 0) = k\ and S'(l, 1) = (/c + 1)! — A;!, which coincides with the initial 
values of the right hand side of f l5.5p . Thus f lS.Sp holds for n > N > n — k. I 

The following identity was found by Schmidt and has been proved in several ways, 
see W^m. 



Example 5.3 (The Apery-Schmidt-Strehl Identity) For n > 0, we have 

n\ f n + k\ fk\^ v-^ f f n + k^ ^ 



,k J \ k I \i I ^ \ k J \ k 



Proof. Let 



n\ I n + k 



and 

n 

= ^fi^^^)9{n, k). 

k=0 

By Zeilberger's algorithm, g{n^ k) is annihilated by the operator 

-8(A; + if - (7A;2 + 2lk + IQ)K + {k + 2fK^. 
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where K denotes the shift operator on k. Now applying the Abel-Zeilberger algorithm 
to /(n, k) and g{n, k), we obtain 

^, , (3 + 2n)(39 + 51n + 17n2)^^ + 2)^ ^ 

S{n) - — -S{n + 1 + 7 7^S{n + 2=0, 

^ ^ (n + 1)3 ^ ^ (n + 1)3 ^ ^ 

Meanwhile, using Zeilberger's algorithm, we find that the Apery numbers 

n\ I n + k 



k 



satisfy the same recurrence relation. Finally, by comparing the initial values, we arrive 
at dsi]). I 

To conclude this paper, we consider the following summations, 

n / k 



k=0 \j=o 



For a = 1,2 and 3, closed forms for S'n"^ have been derived by Andrews and Paule [3| by 
using the method of MacMahon's partition analysis. It is easy to see that these formulas 
can be derived by using the Abel-Zeilberger algorithm. 



Example 5.4 We have 

n k 



k=0 j=o ^-^^ 



(5.7) 



and 



n / k ^2 



k=0 \j=0 



Proof. Let 



It is clear that 



/(n,fc) = l, g{n,k) = J2('') 



g{n,k + l) = g{n,k)+ (^^"^^j and g{n + 1, k) = 2g{n, k) - (^^j. (5.9) 
Applying the Abel-Zeilberger algorithm to f{n, k) and g{n, k), we find that 

5« = E(C l)- T{n + 1) + T(0), 

fc=o V + / 
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where T{k) = (C + 1 — k)g{n, k) and C is a constant. Setting C = —1, we get T(0) = 0. 
Thus 



fc=0 \ "T / 



l)g{n,n + l) 



n / \ n+1 / 

A;=0 ^ ^ j=0 

= -nT-^ + (n + 1)2" 
= n2"-^ + 2". 

Now we consider the evaluation of Sn"* ■ Let f{n,k),g{n,k) be given as above, and 
let h{n, k) = g{n, k)"^. From f l5.9p . we see that 

h{n, k + 1) = h{n, k) + u{n, k), 
h{n + 1, /c) = 4/i(n, /c) + v{n, k), 

where 

.(n,fc)=2^^^J,(n,fc)+^^^^ 

i;(n,A;) = -4 (^^^ ^(n, A;) + {^^ 

It should be mentioned that there is actually no need to impose the condition for u{n^ k) 
and f (n, k) to be hypergeometric in the Abel-Zeilberger algorithm. Therefore, we can 
still apply the Abel-Zeilberger algorithm to /(n, k) and h{n, k) to deduce that 

g{n,k)+ \ I +(r2-C)/i(n,n + l) + (C+l)/i(n,0), 



k=0 \ ^ 



k + 1 



where C is a constant. Setting C = n/2 — 1 and applying the Abel-Zeilberger algorithm, 
we find that 



n /■ \ n 

E(-/2-i-^)U + ib(^'^) = E 

k=0 ^ ^ k=0 



-n + k + 1 



2 



n \ n , 



Hence 



A:=0 ^ ^ k=0 ^ ^ 



^ + 1^2- + ^ 



2 





fc=0 ^ 


n 




^2 


(f: 


/n 






+ l) 2 


V2 





- 1 1 + (| + l") 2^" - ? 



2n ^/2n 



2 V n 
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Using the same argument, we can deduce Calkin's identity [5] 

S^^^ = n23"-i ^ - 3n2""2 (^^^ . (5. 10) 
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